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Abstract 

Probabilities of formation of various hydrogenic species during the exothermic nuclear 
(n, 3 He; t,p)— reaction of atomic helium-3 with slow neutrons are determined. In particular, we 
have found that the probability to form the tritium atom 3 H in its ground state is ~ 55.19287 %, 
while the analogous probability to form the protium atom 1 H in its ground state is ~ 1.02363 %. 
Analogous probabilities of formation of the negatively charged hydrogen ions, i.e. the 3 H" and 1 H~ 
ions, in the nuclear (n, 3 He; t,p) — reaction with slow neutrons, are ~ 7.8680 % and ~ 0.06583 %, 
respectively. We also consider bremsstrahlung from fast fission-type reactions in atomic systems. 
The spectrum of emitted radiation is analyzed. 
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I. INTRODUCTION 

The reaction of the 3 He nuclei with slow neutrons is written in the form [l| 

3 He + n = p + t + 0.764 MeV (1) 

where the notations p and t stand for the protium X H and tritium 3 H nuclei, respectively. 
The cross-section a of this nuclear reaction is extremely large for slow neutrons |2| (a max ~ 
5330- 10~ 24 cm 2 or 5330 barn, for short). The velocities of the two nuclear fragments formed 
in the reaction, Eq.([T]), are w 1.59632 a.u. (tritium) and ~ 4.78897 a.u. (protium). All values 
in this study are given in atomic units, where fi = l,m e = l,e = 1. The unit of atomic 
velocity is v e = ac ~ ~ 2.1882661 ■ 10 8 cm ■ sec' 1 , where c is the speed of light and a is 
the fine structure constant. The atomic velocity v e is the velocity of the Is— electron in the 



hydrogen atom with t 
In our earlier work 



re infinitely heavy nucleus °°H. 

3} we have determined the probabilities to detect various final (bound) 
states in hydrogen atoms which arise during the exothermic nuclear (n, 3 He; t,p)— reactions 
of the one-electron helium-3 ion ( 3 He + ) with slow neutrons. In many actual experiments, 
however, the nuclear (n, 3 He; t,p)— reaction involves the two-electron helium-3 atom rather 
the bare helium-3 nucleus. Analogous nuclear reactions of the 6 Li, 10 B, 14 N nuclei (see, 
e.g., [lj], jsj) with slow neutrons also involve few-electron atoms and ions, rather than the 
bare nuclei and/or one-electron ions. Consequently the products of such reactions may also 
contain a few bound electrons. Our aim below is to evaluate the probabilities of formation 
of such few-electron species. 

In order to determine the probabilities of formation of few-electron atoms and ions during 
nuclear reactions one needs to solve a number of problems. Some of these problems are 
related to the analysis of the basic principles of perturbation theory in quantum mechanics 
4]. In particular, the theory of sudden approximations jlj], (]] plays a central role in our 
analysis below. Based on the sudden approximation we have developed a transparent and 
reliable method which is used to determine the final state probabilities for each possible 
system/state which can arise after the nuclear reaction, Eq.([T]). On the other hand, it is not 
clear what procedure can be used for actual computations of the final state probabilities. In 
[3I we have described a general method based on the use of Nf— electron density matrixes 
for the iVj-electron atomic systems (here Nf < iVj). This is a rigorous method, but in real 
applications it produces a very complex procedure. In many cases this procedure leads to 



the loss of numerical accuracy, e.g., in those cases when one uses highly accurate variational 
wave functions constructed for few-electron atomic systems. Another group of complications 
is related to the electron-electron correlation which cannot be ignored in the incident 3 He 
atom and final H~ ion. 

In this study we have re-considered all mentioned problems and developed an approach 
which can be used in actual calculations of the final probabilities in arbitrary atomic systems 
undergoing nuclear reactions. All details of our analysis will be published elsewhere [7|. Here 
we want to apply this method to obtain the actual data for the nuclear (n, 3 He; t, p)— reaction 
in the two-electron helium-3 atom. A successful test of our method for the 3 He atom with 
the exothermic nuclear (n; t)— reaction, Eq.flTJ), is extremely important. If our method works 
in this case, then it can be applied to other similar reactions in few-electron atomic systems. 
Another goal of this study is to analyze radiation emitted during fast fission-type reactions 
in atomic systems. In particular, we consider the spectrum of bremsstrahlung emitted in 
the nuclear reaction, Eq.flT]), of atomic helium-3 with slow neutrons. 



II. FORMATION OF THE TRITIUM/PROTIUM ATOMS AND IONS 

Note that all current evaluations of the final state probabilities in various atomic systems 

n n 

with nuclear reactions are based on the sudden approximation [5|], |6J. For atomic systems 
with nuclear reactions the sudden approximation j^t means that the electron density of 
the incident atom does not change noticeably during the nuclear reaction in its nucleus. 
In addition to this, the sudden approximation for atomic systems means that the nucleus 
does not move substantially during the nuclear reaction. In general, this is true, if the 
nuclear reaction time r r is significantly shorter than the typical atomic time r a = w 
2.418884331 ■ 10~ 17 sec. It can be shown that these conditions are always obeyed in the case 
of highly exothermic (n; q)— reactions of slow neutrons with light nuclei. This means that 
all changes of electron densities in arising atomic fragments can be described in terms of the 
sudden approximation. 

Our method is also based on the use of the sudden approximation for atomic systems. In 
this work we present only a few basic equations which are used to compute the final state 
probabilities for light atoms and ions. Other details and discussion can be found in [7]. In 
general, for few-electron atomic systems with A— electrons in the incident state we compute 
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the following integral which is called the probability amplitude: 

K 

M in . fi = (* /i (xi,x 2 , . . . ,x#-) | J|exp(zV ■ Ti) | * in (xi,x 2 , . . . ,xjy)) (2) 

i=i 

where is the total number of bound electrons in the final atomic system K < N . Here 
and below the notation Xj = (r», <jj) means the complete set of spatial and spin coordinates 
of the 2-th electron. If the reaction Eq.fJTJ involves the two-electron He atom rather than the 
bare He nucleus, then for K = 1 one finds the 3 H atom (or X H atom) in the final state, while 
K = 2 corresponds to the two-electron 3 H~ ion (or *H ion) in the final state. For K = the 
final system does not contain any bound electrons. The last case (i.e. when K = 0) is not 
of interest for this study. 

In the case of reaction Eq.flTJ in the He atom (where N = 2) with K = 2 one needs to 
compute the following two-electron matrix element (or probability amplitude) 



where V is the velocity of the final nucleus formed in the reaction Eq.([T]). Analogously, for 
the reaction Eq.([T]) in the He atom (where N = 2) with K = 1 we need to determine a 
slightly different two-electron integral (or probability amplitude): 



In Eqs.(j3]) and (j3J) the notation \Ijj n (xi,x 2 ) stands for the two-electron wave function of the 
helium atom, while analogous notations \E f jj(x 1 ,x 2 ) and/or \E^j(xi) designate the two- and 
one-electron wave functions of the final atomic systems, respectively. Below, all two-electron 
wave functions are assumed to be properly antisymmetrized upon all electron variables. 

Let us apply the formulas Eqs.(|3]) and (J4j) to the actual nuclear (n; t)— reaction in the two- 
electron He atom. We shall assume that before such a reaction the He atom was in its ground 
l 1 ^— state. Also, for simplicity, in these calculations all nuclear masses are assumed to be 
infinite. The main interest is to compute the ground-ground state probabilities, since these 
values determine the scale of other similar probabilities. Moreover, if the incident and final 
atomic systems are in their ground states, then the corresponding probability amplitude is 
relatively large. In this case it is very easy to detect and correct possible numerical mistakes. 
The results of our calculations include the corresponding probability amplitudes M in .ji, the 
final state probabilities Pm-ji —\ M in .ji | 2 for both nuclei formed in the reaction Eq. ((TJ) , 



M in .fi = (*/i(xi,x 2 ) | exp(«V -ri + zV- r 2 ) | *m(xi,x 2 )) 



(3) 



M in .ji = (*/<(xi) | exp(zV ■ n) | * in (xi,x 2 )) 



(4) 
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i.e. for the tritium 3 H and protium X H nuclei. Note that the velocities of these two final 
nuclei differ from each other by a factor of ~ 3. This explains very significant differences 
between the final state probabilities determined in calculations for the tritium and protium 
atoms/ions (see below). 

As the starting point for numerical evaluations we can use the following approximate 
expression for the ground l 1 ^— state wave function \l/(xi,X2) of the helium atom (see, e.g., 

B) 

A 3 

*(xi, x 2 ) = — exp(-Ari - Ar 2 ){a(3 - Pa) (5) 

7T 

where the A value is called the effective nuclear charge. For the actual °°He atom one finds 
A = Q — jq, where Q = 2. By performing integration over all spin and angular variables 
one can produce the following formula for the corresponding probability amplitude (for the 
He — > H transition during the nuclear reaction, Eq.flT])) 

M in . fi = 2A~2 / exp(-Ar)j e (Vr)R ne ( ir )r 2 dr (6) 

J 

where je(x) is the Bessel function of the first kind (see, e.g., js], 10]) and Rnii^yr) are the 
radial functions of the hydrogen-like atom/ion (see, e.g., J4]) formed in the final state. If the 
final one-electron atomic system is in its ground state, then we have from Eq.(j6]) 

3 r+°° 

M in , fi = A{A^ ex V [-{A + 1 )r]j G {Vr)r 2 dr (7) 

J 

where V is the velocity of the final fragment, i.e. the velocity of the hydrogen atom. The 
final state probability P gg is 

p -\m i 2 — 64 ^) 3 1 m 
p 99 -I M m;/i 1 - (A + 7)3 f— 7T\¥ (8) 

For the hydrogen atom the parameter 7 = 1, while for the two-electron He atom A — ||. 
In this case for the tritium atom one finds from Eq.(JS]) P gg w 24.37809 %. This value 
corresponds to the process in which one bound a— electron remains bound in the final state. 
Since we have two bound electrons (one a— electron and one (3— electron) in the incident He 
atom, then the final state probability obtained above must be doubled, i.e. we have P gg ~ 
48.74175 %. The use of the highly accurate wave function for the He atom (400 exponential 
basis wave functions [llj) produces for the P gg probability w 55.19287 %. The analogous 
probability for the protium atom : H is ~ 1.02363 %. By using the formulas derived in B 
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and the approximate wave function of the He atom, Eq. (J3J), one can evaluate the final state 
probabilities for some low-lying excited states with n < 3, where n is the principal quantum 
number, in the tritium and protium atoms from the nuclear reaction, Eq.flTJ). These values 
can be found in Table I. 

It is interesting to evaluate the final state probability for the two-electron negatively 
charged hydrogen ion which can be formed during nuclear reaction, Eq.flTJ. To determine 
this probability we used the highly accurate variational wave functions known for the He 
atom and H~ ion 11] . Each of these two functions contains 400 exponents in the relative 



coordinates (basis functions). The total energies produced with such wave functions for the 
ground l 1 ^— states in these two systems are -2.90372437703405 a.u. for the °°He atom and - 
0.527751016544308 a.u. for the °°H _ ion, respectively. By performing numerical calculations 
with these variational wave functions we have evaluated the probability to form the nega- 
tively charged tritium 3 H~ ion in the nuclear (n, 3 He; t,p)— reaction with slow neutrons as ~ 
7.8680 %. The analogous probability to form the negatively charged protium 1 H~ ion was 
evaluated as ~ 0.06583 %. Such a relatively large probability for the newly formed tritium 
ions 3 H~ means that these ions can be detected in modern experiments. Currently, these 
and other similar experiments to detect various atomic species formed during nuclear reac- 
tion Eq.flT]) in two-electron helium-3 atom are critically needed to guide future theoretical 
development. 



III. BREMSSTRAHLUNG FROM FISSION-TYPE PROCESS IN ATOMIC SYS- 
TEMS 

Another group of problems is related to the analysis of radiation emitted during the 
reaction of the 3 He nuclei with slow neutrons, Eq.flTJ). It is clear that such radiation must be 
similar to radiation emitted during an arbitrary fission-type reaction in few-electron atomic 
systems (atoms, ions, etc). First, let us consider a fast nuclear fission-type reaction in a 
one-electron atomic system. During a fission-type nuclear reaction the electron becomes 
free, and will interact with the rapidly moving nuclear fragments. Such an interaction will 
produce radiation which is, in fact, a breaking radiation or bremsstrahlung. Here we want 
to derive the explicit formulas for this radiation and investigate its spectrum. The electric 
charges and velocities of these fission fragments are Qie,Q 2 e and Vi,V 2 , respectively. To 
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simplify all formulas below, we shall assume that the both fission fragments move along the 
Z— axis. Furthermore, in this study we restrict ourselves to the consideration of the non- 
relativistic processes only. Briefly, this means that the two velocities V\, V 2 are assumed to 
be significantly smaller than the speed of light c. The case of arbitrary velocities is discussed 
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12|. 



The second time-derivative of the dipole moment d is d = er, where r = (x, y, z) is the 
radius-vector of the accelerated electron. In the case of a fission-type reaction in atomic 
systems the explicit formula for electron's acceleration is 



1 r 

rn 



Q 2 e\i Qie 



1 

^R 2 



(9) 



where R\ = J x 2 + y 2 + (z — Vit) 2 and R 2 = \Jx 2 + y 2 + (z + V 2 t) 2 . The notations V\ and V 2 
stand for the velocities of the first and second fission fragments. The second time-derivative 
of the dipole moment is 



9^ v( ±) + 9*l v <± 



Q ie 2 Rx Q 2 e 2 R 2 



m e Rf 



m P Rn 



(10) 



where Ri = (x, y,z — Vit) and R 2 = (x, y,z + V 2 t) are the three-dimensional vectors. The 
intensity of the non-relativistic bremsstrahlung from a fission-type reaction in a one-electron 
atomic system is 



dl = ^—(dxn) 2 dn= (—) 



Atcc 3 



e 



w^g"^Rf ~q ~rJ 



sin 9dQ 



(11) 



where 9 is the angle between the vector d and the vector n which designates the direction of 
propagation of radiation. The notation Q in the last equation stands for an arbitrary electric 
charge. This value can be considered as an additional input parameter. In particular, one 
can choose Q = Q\, or Q = Q 2 . In the case of the fast nuclear fission- type reaction in a 
few-electron atomic system the last formula takes the form 



dl 



e 2 \2N e Q 2 Ri Q2 R2 
Q Rf Q Ri 



(-)' 



Atcc 3 



sin 2 6dVt 



(12) 



where N e is the number of electrons which become free after this nuclear reaction. The 
formula, Eq. f|T2|) . can be re- written in the form 

, „2 . „ at r>2 _ra 1 ra 



dQ 



e 2 ^2N e Q 2 { QiJ_ QiJ_ giga/Ri-Ra 
Q 2 Rf Q 2 RI Q 2 



(-) 



47TC 3 



■(- 

[ RlRl 



sin 9 



(13) 
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for the differential cross-section. Our derivation of these formulas for few-electron systems is 
based on the fact that radiation emitted by different post-atomic electrons is non-coherent. 
Indeed, different free electrons move as independent (or 'random') quantum particles. Note 
that the bremsstrahlung from fast fission-type process in atomic systems can be coherent 
at certain conditions, e.g., if two fission fragments move with very large velocities, while all 
post-atomic electrons almost do not move at the beginning of the process. In this case, one 
needs to introduce an additional factor N e in the last formula, Eq. fflU]) . 

Let us discuss the spectrum of the emitted radiation. As follows from the formula, 
Eq.f lTSj) . the intensity of bremsstrahlung from fission-type processes in atomic systems 
rapidly decreases with the time I ~ t~ 4 after the nuclear reaction which proceeds at t = 0. 
Such a behaviour is typical for atomic systems with fast nuclear reactions. The spectral 
resolution R(u) (or spectral function) of the intensity of dipole radiation is written in the 
form 

s 4 •• , 9 duo 4a; 4 . , ,9 dcj , JS 

^ = 3? I *• I' S = 3? I I S (14) 
where is the Fourier component of the dipole moment d introduced above. As follows 
from Eq.( ll4p . to determine the spectral resolution R(u) of bremsstrahlung from a fission- 
type process one needs to find the Fourier components of the dipole moment. In this study 
we apply a different method and determine the Fourier components of the d w vector, Eq.Q. 
Finally, the problem is reduced to the calculation of the two following Fourier transformations 

I\ (uj ! ~j V, cos tj) = / (— ± 2— cost] ■ t + t J 2 exp(iut)dt (15) 
V Jo K V Z V ' 

and 

d r+oo . q2 q 3 Q q 

h(u;y,V,cosr)) = J y— ± 2— cosrj ■ t + t J Hexp(iut)dt = -i—I^u] —,V,cosr]) 

(16) 

where a is the 'effective' radius of the original electron shell and rj is the angle between 
the electron acceleration and Z— axis (or the line of the nuclear motion). Note that the 
spectrum of the emitted radiation depends explicitly upon cos?]. This dependence can be 
found for all fast fission-type reactions and processes. The lower limit in these formulas is 
zero, but, in reality, we cannot use times t which are shorter than r = -. Formally, this 
means that the lower limits in Eqs. f|T5|) and f lT6|) can slightly be changed and this can be 
used to simplify the explicit formulas for the corresponding Fourier components. The upper 
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limits in Eqs. (ll5[) and ( IT6|) can also be chosen to be finite. For instance, it is possible to 
obtain a very good approximation for the integrals, Eqs.f lT5|) and f|T6|) . by using the upper 
limit In actual situations one can use the formula exp(zx) = cosx + isinx and then 
calculate all arising integrals numerically, or by using analytical formulas from the Tables of 
Fourier sine /cosine transformations (see, e.g., 13[). A good analytical approximation for 
the integral, Eq.( !T5|) . is 

77 , V, cos?]) ~ exp(— no — cos m / [b 2 + t 2 ]~2 exp(zur)dT 
V ^ V 'Jo 

I a \ yuj T ^ . 1 A + ^bu\ 7r u 2 i /irA 

= exp — uo — cos n — KAbu) H — in f- — ^ (17) 

where 6 = anc [ ^(x) is the MacDonald function (see, e.g., The derivative of 

Eq. ffTTl) with respect to the frequency u allows one to determine the ^(^j "7? ^ cos 77) integral, 
Eq.( !T6|) . and obtain the analytical formula for the spectral resolution R(u) of the intensity of 
bremsstrahlung, Eq. ([T4]) . from the fast fission-type (fast) process in few- and many-electron 
atomic systems. The derivation of explicit expression for the spectral resolution R(u) is 
straightforward, but the final formula is cumbersome and it is not given here. Note only 
that the term ~ Q1Q2 in this formula represents interference of radiation which is emitted 
by the two rapidly moving fission fragments. In application to the nuclear fission of heavy 
nuclei, e.g., 239 Pu, 247 Cm, etc, the spectral resolution R(u) must be averaged over the actual 
mass/charge distribution known for each fissionable nucleus. 



IV. CONCLUSION 

We have considered the nuclear reaction, Eq.([T]), involving the two-electron 3 He atom 
and slow neutrons. To determine the final state probabilities for one- and two-electron 
tritium and protium atoms/ions we have developed a new method which is based on direct 
computation of probability amplitudes. This method is much simpler than an alternative 
procedure based on the use of few-electron density matrices. Furthermore, it allows one 
to determine the final state probabilities for various bound states in the incident and final 
atomic systems. By using this method we have found that the probability to form the 
tritium atom 3 H in its ground state is ~ 55.19287 %, while analogous probability to form 
the protium atom in its ground state 1 H is ~ 1.02363 %. Analogous probabilities for the 
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negatively charged hydrogen ions, i.e. for the 3 H~ and 1 H~ ions to be formed in the nuclear 
(n, 3 He; t,p)— reaction with slow neutrons, are ~ 7.8680 % and ~ 0.06583 %, respectively. 
The corresponding uncertainties for such probabilities can be evaluated as ~ 0.3 — 1 • 10 -4 
% for the ground states and as ~ 0.5 — 1 • 10 -3 % for the excited states. Another important 
problem is the analysis of bremsstrahlung from fast fission-type reactions in atomic systems. 
In this study such an analysis is performed for the non-relativistic case. Its generalization 
to the systems with relativistic velocities is quite complicated (see, e.g., [l^j]). 

In the future it is necessary to study analogous nuclear reactions involving atoms with 
the 6 Li, 10 B, 14 N (and some other) nuclei. For such atomic systems we need to determine the 
final state probabilities to form various few-electron atoms and ions in the reaction of these 
atoms with slow neutrons. These results are of interest in numerous applications of these 
nuclear reactions, e.g., in the Boron Neutron Capture Therapy (BNCT) where the reaction 



[n, 10 B; 7 Li, 4 He) is used [15| - [17|. Nuclear reaction Eq.([T]) and the analogous reaction with 
the 6 Li nuclei are extensively used in thermonuclear explosive devices [3] ■ Also, these nuclear 
reactions are used in various detectors (see, e.g., [l^) of thermal and slow neutrons. 

By using our results for the final state probabilities of different atomic species formed in 
such reactions we can improve the overall sensitivity of such detectors. 
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TABLE I: The probabilities (in %) of the final states in the tritium/protium atom arising in the 
exothermic nuclear reaction, Eq.(l), of the two-electron helium-3 atom with slow neutrons. Here 
the notation a(b) stands for a ■ 1CT 6 . 

atom/state Is 2s 2p 3s 3p 3d 

3 H 48.7418 9.42077 0.803701 0.534751 0.053371 0.032033 

l R 0.53715 0.105205 0.002480 0.008709 1.28725-10" 2 3.7800-1Q- 6 
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